function [] = fig_sec11_4_congControlAndBackpressureRouting_dualDecomp()

close all;

hd = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_hd.txt');
xp = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_xp.txt');
qnd = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_qnd.txt');
objFunc = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_objFunc.txt');

jom_hd = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_jom_hd.txt');
jom_xp = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_jom_xp.txt');
jom_qnd = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_jom_qnd.txt');
jom_objFunc = load ('data/crossLayerCongControlBackpressureRoutingDualDecomp_jom_objFunc.txt');

multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-2 0] , '-' , 'Q_{nd}' , qnd , [] , [0 1300] , '-' , 'x_p' , xp , jom_xp , [0 26] , '-' , 'h_d' , hd , jom_hd , [0 26] , '-');


jom_objFunc
grad_objFunc = objFunc (end)

jom_qnd
qnd (end,2:end)
